Influence of Roughness and Disorder on Tunneling Magnetoresistance 
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A systematic, quantitative study of the effect of interface roughness and disorder on the magne- 
toresistance of FeCo|vacuum|FeCo magnetic tunnel junctions is presented based upon parameter-free 
electronic structure calculations. Surface roughness is found to have a very strong effect on the spin- 
polarized transport while that of disorder in the leads (leads consisting of a substitutional alloy) is 
weaker but still sufficient to suppress the huge tunneling magneto-resistance (TMR) predicted for 
ideal systems. 



Tunneling magnetoresistance (TMR) refers to 
the dependence of the resistance of a FMi|I|FMa 
(ferromagnet | insulator | ferromagnet) magnetic tunnel 
junction (MTJ) on the relative orientation of the magne- 
tization directions of the ferromagnetic electrodes when 
these are changed from being antiparallel (AP) to paral- 
lel (P): TMR = {Rap - R P )/R P = (G P - G A p)/G A p- 
Since the discovery of large values of TMR in MTJs 
based upon ultrathin layers of amorphous AI2O3 as 
insulator^ a considerable effort has been devoted to 
exploiting the effect in sensors and as the basis for 
non-volatile memory elements. Understanding TMR 
has been complicated by the difficulty of experimentally 
characterizing FM|I interfaces. The chemical composi- 
tion of the interface has been shown£ to have a strong 
influence on the magnitude and polarization of the TMR 
and knowledge of the interface structure is a necessary 
preliminary to analyzing MTJs theoretically. In the 
absence of detailed structural models of the junctions 
and the materials-specific electronic structures which 
could be calculated with such models, the effect was in- 
terpreted in terms of electrode conduction-electron spin 
polarizations Pi, using a model suggested by Julliere 3 
in which the TMR = 2P 1 P 2 /(1 - PiP 2 )- A great deal 
of discussion has focussed on the factors contributing 
to the quantitj«* P but the use of amorphous oxide as 
barrier material made impossible a detailed theoretical 
study with which to confront experiment 

The situation changed quite drastically with the recent 
observation of large values of TMR at room temperature 
in FeCo|MgO|FeCo MTJs in which the MgO tunnel bar- 
rier was mono-** or poly-crystalline,' 1 - This work was 
motivated in part by the predictionii^ by materials- 
specific transport calculations of huge TMR values for 
ideal Fe|MgO|Fe structures. This new development lends 
fresh urgency to the need to understand the factors gov- 
erning the sign and magnitude of TMR because the 
largest observed value of 353% at low temperature^ is 
still well below the ab-initio predicted values of order 
10,000% for the relevant thicknesses of MgO^i Some ef- 
fort has been devoted to explaining the discrepancy in 
terms of interface relaxation^ 3 " or the formation of a layer 



of FeO at the interfac e - 1 15 but the role of interface dis- 
order has only been speculated upon^& 

Method . In this paper, we use first principles elec- 
tronic structure calculations to study the effect of rough- 
ness and alloy disorder on TMR in MTJs with a vac- 
uum barrier and Fei-^Coa; alloy electrodes. A vac- 
uum barrier was chosen for its simplicity and because 
there are many studies of spin-dependent vacuum tun- 
neling in its own right. 17,1 -8il2i2Bi21 We consider the ef- 
fect of diffusive scattering in the linear-response regime 
in a two-step procedure. In the first step, the elec- 
tronic structure of the Fei-zCozlvacuumlFei-zCoz MTJ 
is determined using the local-density approximation 2 ^ 
of density functional theory. The self-consistent cal- 
culations are performed with the tight-binding linear 
muffin-tin orbital (TB-LMTO)^ 3 . surface Green's func- 
tion method" and disordered systems are treated using 
the layer CPA (coherent potential approximation) 2^ The 
atomic sphere (AS) potentials serve as input to the sec- 
ond step in which the transmission matrix entering Lan- 
dauer's transport formalism 2 ^ is calculated using a TB- 
MTO implementation 2 !^ of a wave-function matching 
scheme due to Ando. 29 Disorder is modelled in large lat- 
eral supercells by distributing the self-consistently calcu- 
lated CPA-AS potentials at random, layer-for-layer in the 
appropriate concentrations for as many configurations as 
are required. In most of the calculations to be presented, 
supercells containing 10 x 10 atoms per monatomic layer 
were used. 

We consider transport in the (001) growth direction 
keeping atoms at the surfaces unrelaxed in their bulk bec 
positions. For Fe leads the experimental lattice constant 
ape = 2.866A is used. For Co leads oc = 2.817A is 
chosen so that the bec lattice has the same volume as 
hep Co. The alloy is considered to obey Vegard's law 
whereby 

aF Cl _ x Co x = (1 - x)a Fc + xa Co (1) 

The vacuum region is modelled in the atomic spheres ap- 
proximation (ASA) by filling the space between the leads 
with 'empty' spheres of the same size 3 ^ and on the same 
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FIG. 1: (Color online) Conductances Gp in (▼), Gp aj (A), 



and G° 



of an Fel vacuum IFe MTJ as a function of the 



barrier thickness (measured in units of layers of a bec lattice). 
The solid lines are for ideal junctions. The dashed lines are 
configuration-averaged conductances for rough junctions pre- 
pared by removing, at random, half of the Fe atoms from one 
surface and depositing them, at random, on the other sur- 
face. The results are normalized to the lxl surface unit cell 
used for the ideal case. The large symbols refer to data points 
which also appear in the next figure. 



bec lattice as in the leads. Consequently we measure the 
thickness of the vacuum barrier in monolayers (MLs) of 
empty atomic spheres. 

Ideal Fe|vac|Fe. The conductances G P and G a AP for 
an ideal, ordered Fe| vacuum |Fe MTJ are shown in Fig. 2] 
(solid lines) as a function of the width of the barrier for 
the minority and majority spin channels, a = min,maj. 
In all three cases an exponential dependence on the bar- 
rier width is reached asymptotically. For Gp aj , the 
asymptotic dependence sets in after about five MLs of 
vacuum; the initial subexponential behaviour is related 
to the deviation of the barrier from a simple step form. 
The behaviour of Gp m is much more complex, becoming 
exponential only when the vacuum is some 8 MLs or more 
thick; the AP case is intermediate. The TMR increases 31 
as a function of the width of the vacuum barrier reaching 
huge asymptotic values of order 20,000%, similar in size 
to those reporteciii*i£ for Fe|MgO|Fe. The non-trivial de- 
pendence of Gp m on barrier width can be understood in 



terms of a surface stated at T (k|| = 0) which is very 
close to the Fermi level in the minority spin channel of 
iron but well below it for majority spin electrons. In 
the P configuration one such state exists on each sur- 
face and at values of k|| ^ these states become surface 
resonances which form bonding-antibonding pairs with 
transmission probabilities close to unity^ As the vac- 
uum barrier is widened, the coupling between the surface 
resonances weakens until the bonding-antibonding split- 
ting becomes smaller than the resonance width at which 



point the maximum transmission becomes smaller than 
unity and normal exponential dependence of the conduc- 
tance on the barrier width sets in; the asymptotic slope in 
Fig. n is consistent with our calculated Fe workfunction 
of 4.7 eV. 

Roughness: 50% coverage. An ideal tunnel junction 
such as that just considered is impossible to realize in 
practice; there will always be some finite amount of disor- 
der whether it be surface roughness, islands, dislocations 
etc. A more realistic model is obtained by considering 
an Fe|vacuum|Fe system with incomplete (rough) surface 
layers, modelled by occupying, at random, a fraction of 
the lattice sites of the topmost layer with Fe atoms (see 
the inset in Fig. 1 ) and using the layer CPA to determine 
the corresponding AS potentials. 

Roughness resulting from depositing some Fe atoms on 
an ideal surface has two effects. The first is to destroy 
the point group symmetry which led to the existence of 
a symmetry gap at T for states with Ai symmetry, and 
thus the surface state with that symmetry. We expect 
this to reduce the large minority spin conductance found 
for the ideal vacuum barrier. The second effect of rough- 
ness is to reduce the width of the vacuum barrier and thus 
to enhance the conductance which depends exponentially 
on the barrier width. 

These competing effects are disentangled by starting 
with an ideal vacuum barrier and moving half of the Fe 
atoms from one surface and depositing them on the other 
so that the barrier width is, on average, unchanged. The 
results of these supercell calculations, averaged over 20 
configurations and normalized to the ideal MTJ (lxl 
surface unit cell) results are included in Fig.^ The con- 
tribution of the interface states to G™ m is completely 
quenched by the roughness and the conductance reduced 
by four orders of magnitude. Up is enhanced by about 
an order of magnitude because the increased conductance 
from those parts of the rough surfaces which are closer 
than average more than compensates the decreased con- 
ductance from part of the rough surfaces which are fur- 
ther away than average. All conductances now exhibit 
the same qualitative behavior: exponential decay. Not 
only is the absolute value of the TMR greatly reduced 
by roughness, the sign of the polarization is reversed. 

This result immediately prompts us to ask how the 
TMR depends on the amount of roughness and, in par- 
ticular, what coverage is needed to suppress the contri- 
bution of the resonances to Gp m . This issue is addressed 
in Fig. |21 where the conductance of a MTJ is shown as a 
function of surface coverage. Zero coverage corresponds 
to an ideal MTJ with 8 MLs vacuum; 100% coverage (not 
shown) to 6 MLs. We see that 5% coverage is sufficient to 
reduce G P im , and the TMR, by two orders of magnitude. 
Unless the surface roughness is less than a few percent, 
the TMR lies between about 10 and 1000% (Fig.0 inset), 
comparable to values observed in experiment. As G™ m 
is reduced by roughness, Gp 10 " 3 increases monotonically 
with increasing coverage as the average barrier width de- 
creases; the two cross at a surface coverage of about 25%. 
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FIG. 2: (Color online) Configuration-averaged conductances 
Gp in (T), Gp aj (A), and G AP (x) of an Fe|vacuum|Fe MTJ 
with 8 ML thick vacuum barrier as a function of the surface 
coverage, normalized to a 1 x 1 surface unit cell. The dashed 
line denotes G AP predicted from Eq. [5] Large symbols refer 
to the similarly marked data points in the previous figure. 
Inset: TMR as a function of the surface coverage. The dashed 
fine is the value predicted using Julliere's expression and a 
calculated DOS polarization of 55%. 



FIG. 3: (Color online) G^ ln (T), Gp aj (A), and G AP (x) 
for an Fei-^Co^ IvacuumlFei-xCoz MTJ with 8 ML thick 
vacuum barrier as a function of x, the concentration of Co 
atoms, calculated in the virtual crystal (VCA: solid lines) and 
CPA/supercell (SC: dashed lines) approximations. Conduc- 
tances are configuration averaged and normalized to a 1 x 1 
surface unit cell. Inset: TMR in VCA (o) and CPA/SC (•) 
approximations. 



The single spin AP conductance, G AP , is described 
qualitatively very well for all coverages by the relation^ 

G AP = ^Gp a3 G™ m (2) 

(dashed line in Fig. [2J and quantitatively for coverages 
greater than a few percent i.e., as soon as the surface reso- 
nance contribution is killed by roughness^ When Eq. (0 
holds, G P = Gp aj + Gp™ is greater than Gap = 2G AP 
and the TMR is always positive although the polariza- 
tion in the P configuration changes sign, reaching its 
minimum value at 25% coverage where Gp^ 3 and Gp m 
crossover. 

As long as there is enough roughness to quench the 
surface resonance, there is order-of-magnitude agreement 
between our TMR calculated as a function of surface 
roughness and the value obtained using Julliere's expres- 
sion for the TMR with the calculated, bulk density-of- 
states polarization for Fe, P = 55%, shown as a dashed 
line in the inset to Fig. [21 

Disorder in the leads . It is interesting to study an 
intermediate type of disorder where there is no rough- 
ness but the electrodes are made of a substitutional 
7e\- x Co x magnetic alloy. As the alloy concentration 
x is increased, the disorder increases but the underly- 
ing electronic structure is also changing as the Fermi 
energy rises. To distinguish these two effects, we first 
carry out calculations for a VCA | vac | VCA MTJ in which 
the Fei-zCoz electrodes are treated within the virtual 



crystal approximation (VCA). This allows us to probe 
the effect of changing the Fermi energy without includ- 
ing disorder. It turns out to be also very convenient 
to embed Fei-^Co^lvaclFei-^Co^ between VCA leads 
(schematically denoted) as VCA|CPA|vac|CPA|VCA for 
performing the self-consistent potential calculations and 
as VCA|SC|vac|SC|VCA for the conductance calculation 
where, as before, CPA potentials are used as input to 
the supercell (SC) transport calculations. Details of the 
thicknesses of the various regions needed to achieve con- 
verged results as well as the number of k-points, supercell 
sizes and other technical details will be given in a forth- 
coming publication. 

The results of these calculations are shown in Fig. [3] 
Gp 10,3 and G AP are unchanged by disorder within the ac- 
curacy of the calculation. The largest change can once 
again be seen in Gp m . In the case of ideal VCA elec- 
trodes, there are localized surface states in the minority 
channel as one would expect. With increasing Co concen- 
tration x, the Fermi level rises with respect to these sur- 
face states, resulting in a decreasing contribution to the 
conductance from the surface resonances. This is clearly 
demonstrated by the behavior of Gp m in Fig. [21 Lead 
disorder quenches the contribution of the resonant states 
by destroying the point group symmetry, eliminating the 
symmetry gap and broadening the resonances. At a con- 
centration of about 40%, the surface resonances are no 
longer dominant and the effect of disorder is quite small; 
for x > 0.4 the trend as a function of x is described quite 
well by the VCA calculation. Eq. © is again found to 
describe the AP conductance very well. We see thus that 
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lead disorder has the same effect as roughness in quench- 
ing the huge values of TMR (inset to Fig. |3J) found for 
the disorder-free VCA reference system, albeit less effec- 
tively. 

Summary. Using first-principles calculations we find 
that both roughness and alloy disorder quench the TMR 
for ideal magnetic tunnel junctions to values compara- 
ble to the highest found experimentally for monocrys- 
talline barrier materials or vacuum tunneling. We there- 
fore propose that these experiments are still in the rough- 
ness/disorder limited regime. 

This work is part of the research program of the 
"Stichting voor Fundamentcel Onderzoek der Materie" 



(FOM) and the use of supercomputer facilities was spon- 
sored by the "Stichting Nationale Computer Facilitcitcn" 
(NCF), both financially supported by the "Nederlandse 
Organisatie voor Wetenschappelijk Onderzoek" (NWO). 
VMK is supported by "NanoNed", a Nanotechnology 
programme of the Dutch Ministry of Economic Affairs. 
KX acknowledges support from NSF of China, Grant 
No. 90303014. 

We are grateful to Anton Starikov for permission to 
use his version of the TB-MTO code based upon sparse 
matrix techniques with which we performed some of the 
calculations. 



1 J. S. Moodera, L. R. Kinder, T. M. Wong, and R. Meser- 
vey, Phys. Rev. Lett. 74, 3273 (1995). 

2 J. M. De Teresa, A. Barthelemy, A. Fert, J. P. Contour, 
F. Montaigne, and P. Seneor, Science 286, 507 (1999). 

3 M. Julliere, Phys. Lett. A 54, 225 (1975). 

4 Pi was not defined strictly in Ref. y and it has frequently 
been interpreted in terms of the density-of-states polar- 
ization P = [Df (ef) - Di(e/)]/[D T (E/) + D L (s f )] where 
D a {e) is a bulk density of states. It has been pointed ouii- 
that other measures of the spin polarization may be more 
appropriate when studying transport. 

5 I. I. Mazin, Phys. Rev. Lett. 83, 1427 (1999). 

6 E. Y. Tsymbal, O. N. Mryasov, and P. R. LeClair, J. Phys.: 
Condens. Matter. 15, R109 (2003). 

7 X. G. Zhang and W. H. Butler, J. Phys.: Condens. Matter. 
15, R1603 (2003). 

8 S. Yuasa, T. Nagahama, A. Fukushima, Y. Suzuki, and 
K. Ando, Nature Materials 3, 868 (2004). 

9 S. Yuasa, T. Katayama, T. Nagahama, A. Fukushima, 
H. Kubota, Y. Suzuki, and K. Ando, Appl. Phys. Lett. 
87, 222508 (2005). 

10 S. S. P. Parkin, C. Kaiser, A. Panchula, P. M. Rice, 
B. Hughes, M. Samant, and S. H. Yang, Nature Materi- 
als 3, 862 (2004). 

11 W. H. Butler, X.-G. Zhang, T. C. Schulthess, and J. M. 
MacLaren, Phys. Rev. B 63, 054416 (2001). 

12 J. Mathon and A. Umerski, Phys. Rev. B 63, 220403 
(2001). 

13 D. Wortmann, G. Bihlmayer, and S. Bfiigel, J. Phys.: Con- 
dens. Matter. 16, s5819 (2004). 

14 X. G. Zhang, W. H. Butler, and A. Bandyopadhyay, Phys. 
Rev. B 68, 092402 (2003). 

15 C. Tusche, H. L. Meyerheim, N. Jedrecy, G. Renaud, 
A. Ernst, J. Henk, P. Bruno, and J. Kirschner, Phys. Rev. 
Lett. 95, 176101 (2005). 

16 K. D. Belashchenko, J. Velev, and E. Y. Tsymbal, Phys. 
Rev. B 72, 140404 (2005). 

17 J. A. Stroscio, D. T. Pierce, A. Davies, R. J. Celotta, and 
M. Weinert, Phys. Rev. Lett. 75, 2960 (1995). 

18 S. F. Alvarado, Phys. Rev. Lett. 75, 513 (1995). 

19 S. N. Okuno, T. Kishi, and K. Tanaka, Phys. Rev. Lett. 



88, 066803 (2002). 

20 H. F. Ding, W. Wulfhekel, J. Henk, P. Bruno, and 
J. Kirschner, Phys. Rev. Lett. 90, 116603 (2003). 

21 M. M. J. Bischoff, T. K. Yamada, C. M. Fang, R. A. de 
Groot, and H. van Kempen, Phys. Rev. B 68, 045422 
(2003). 

22 U. von Barth and L. Hedin, J. Phys. C: Sol. State Phys. 
5, 1629 (1972). 

23 O. K. Andersen, Z. Pawlowska, and O. Jepsen, Phys. Rev. 
B 34, 5253 (1986). 

24 I. Turek, V. Drchal, J. Kudrnovsky, M. Sob, and P. Wein- 
berger, Electronic Structure of Disordered Alloys, Surfaces 
and Interfaces (Kluwer, Boston-London-Dordrecht, 1997). 

25 P. Soven, Phys. Rev. 156, 809 (1967). 

26 S. Datta, Electronic Transport in Mesoscopic Systems 
(Cambridge University Press, Cambridge, 1995). 

27 K. Xia, P. J. Kelly, G. E. W. Bauer, I. Turek, J. Ku- 
drnovsky, and V. Drchal, Phys. Rev. B 63, 064407 (2001). 

28 K. Xia, M. Zwierzycki, M. Talanana, P. J. Kelly, and 
G. E. W. Bauer, Phys. Rev. B 73, 064420 (2006). 

29 T. Ando, Phys. Rev. B 44, 8017 (1991). 

30 The Fe and Co lattice constants correspond to Wigner- 
Seitz atomic sphere radii of 2.667 and 2.621 Bohr atomic 
units, respectively. 

31 The conductance polarization changes sign from positive 
to negative as the vacuum separation is increased. The 
positive polarization found by MacLaren et alw— can be 
reproduced by using non selfconsistent potentials and small 
number of k points. 

32 J. M. MacLaren, X.-G. Zhang, and W. H. Butler, Phys. 
Rev. B 56, 11827 (1997). 

33 O. Wunnicke, N. Papanikolaou, R. Zeller, P. H. Dederichs, 
V. Drchal, and J. Kudrnovsky, Phys. Rev. B 65, 064425 
(2002). 

34 K. D. Belashchenko, E. Y. Tsymbal, M. van Schilfgaarde, 
D. A. Stewart, I. I. Oleynik, and S. S. Jaswal, Phys. Rev. 
B 69, 174408 (2004). 

35 For ideal MTJs, the resonant tunneling contribution will 
be quenched for sufficiently thick barriers so Eq. can 
also be expected to hold in that situation. 



